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ABSTRACT 

We present the results of optical and infrared photometry of the quiescent X-ray tran- 
sient V404 Cyg during the period 1992-2003. The ellipsoidal modulations extracted 
from the most complete databases (years 1992, 1998 and 2001) show unequal maxima 
and minima with relative strength varying from year to year although their peak to 
peak amplitudes remain roughly constant at 0.24±0.01 magnitudes. Fast optical vari- 
ations superimposed on the secondary star's double-humped ellipsoidal modulation 
were detected every year with a mean amplitude of ~0. 07 mags. We have not found 
significant changes in the activity during this decade which indicates that this vari- 
ability is probably not connected to the 1989 outburst. We have found periodicities 
in the 1998 and 2001 data near the 6hr quasi-periodicity observed in 1992, although 
we interpret it as consequence of the appearance of a flare event almost every night. 
Significant variability is also present in the / and near infrared (J and K s ) bands 
and this decreases slightly or remains approximately constant at longer wavelengths. 
A cross correlation analysis shows that both the R and / emission are simultaneous 
down to 40 s. 
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1 INTRODUCTION 

Low-Mass X-ray Binaries (LMXBs), are semi-detached 
binary systems where a normal low-mass late type star 
transfers gas onto a compact remnant i.e. a neutron star 
or black hole. A significant fraction of LMXBs are soft 
X-ray transients (SXTs) characterized by episodic X-ray 
outbursts where the luminosity increases dramatically. 
During the long intervals (usually many years) between 
outbursts, the SXTs are in quiescence and the optical 
emission is dominated by the radiation from the companion 
star. This offers a superb opportunity to analyze and 
obtain dynamical information which eventually enables 
us to constra i n the nature of the compact object (see eg. 
ICharleslll99i ICasares et al.l l200lL The majority of these 
systems appear to contain black holes and they are, in 
fact, the best evidence we have for the existence of stellar 
mass black holes. The most secure case is the SXT V404 
Cyg which has a mass function of 6.08 ± 0.06 Mq (Casares, 
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Charles & Naylor 1992), twice t he usually accepted uppe r 
limit to that of a neutron star llRhoades fc Ruffinill97lh . 
The estimated mass for the compact object is 12±2 M0 



Shahbaz et aTlll994h which, together with GRO 1915+105 
Greiner. Cubv fc McCaughrearJliuOlh . are the most mas- 



sive stellar mass black holes yet measured. 



Furthermore, V404 Cyg is considered a peculiar SXT 
because it has on e the longest orbital periods (6.5 d; 
ICasares et al] 0.992) and is the most X-ray luminous i n 
quiescence (L x ~ 10 33 - 10 34 ergs" 1 : ICarcia et ail 1200 ill . 
Apart from this, it is largely known that V404 Cyg exhibits 
short-term (sub-orbital ) variations during quiescence in 
the o ptical cont inuirm_J^^gner^^l|h^92£ Pavlenko^^l 
| 1996|) and Ha flux JCasares fc CharlesUl992l 1 Casar es et al 
1993) wher e the Ha profiles s eems to be cor related to the 
conti nuum llHvnes e t al]l20o 3l. The infrared llSanw al et al] 
Il995l) . X-ravs iWagner et al.lll99l iKong et al.ll2002ll and 
radio iHiellmiriT^tHlrtTH^OOOf) are also variable. Recent 
simultaneous observations in X-rays and optical show that 
optical line and continuum variability are well correlated 
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with X-rays flares (Hynes et al. in preparation). On the 
other hand it seems that fast variability during quiescence 
is not a peculiarity of V404 Cyg but a fingerprint of 
quiesce nt SXTs, irrespective of the nature of the compact 
object JZurita. Casares fc Shahba3l2003t iHvnes et alJl200l 
hereafter ZCS). Although these variations seem to originate 
in the accretion disc ra ther than in the op tical c ompanion or 
gas s tream/hot spot iZurita fe Casareall998t IHvnes et alJ 
120021 . ZCS) the underlaying physical mechanism is not 
yet understood. The dynamical (Keplerian) timescales 
associated with these fluctuations suggest that the flares 
are produced in the outer disc (ZCS). Alternatively , based 
on the similarities of the power density spectr um (PDS) of 
A0620 -00 to those of low/hard state sources IHvnes et alJ 
(2003) have pointed out that the optical variability is at 
least partly related to the inner accretion flow even in 
quiescence. This suggestion was also supported by the 
discovery of the correlated optical and X-ray flares (Hynes 
et al. in preparation). 

Since V404 Cyg is the brightest SXT, and hence acces- 
sible with 1-m type telescopes, we embarked on a study of 
the fast variability over the period 1992-2003. This paper 
comprises of a simultaneous optical and infrared photomet- 
ric study, with the aim to determining the colour of the 
variability and its possible origin. 



2 OBSERVATIONS AND DATA REDUCTION 

We obtained CCD photometry of V404 Cyg with various 
1 m-class telescopes during the years 1992, 1998, 1999, 
2001, 2002 and 2003. The 1992 data were obtained in 
the Johnson R band with the GEC CCD camera on 
the 1 m JKT telescope at Observatorio del Roque de los 
Muchachos. V404 Cyg was observed in 1998 simultane- 
ously in the R and / band with identical Thomson CCD 
cameras at both the 0.8 m IAC80 and lm OGS telescopes 
at Observatorio del Teide. We also obtained R band 
photometry in 1999 with the SITe2 CCD camera on the 
JKT. Further de tails on the 199 2 and 1999 observations 
can b e found in IPavlenko et~aH l)l99rl) and IHvnes et alJ 
(2002) respectively. During Jul-Sep 2001 -R-band and J and 
K a simultaneous observations were performed using the 
Thomson CCD camera at the 0.8 m IAC80 and the infrared 
camera at the 1.5 m Telescopio Carlos Sanchez (TCS) also 
at Observatorio del Teide. Finally we observed V404 Cyg 
during May-Jun 2002 and Jul 2003 also in the .R-band 
and with identical Thomson CCD cameras at both the 
0.8m IAC80 and lm OGS telescopes at Observatorio del 
Teide. Full details of these observations are given in Table 

All images were corrected for bias and flat-fielded in the 
standard way using iraf. Seeing conditions were not good 
enough to cleanly separate the contribut ion of the target and 
its ne arby (1.4arcsec) line-of-sight star llUdalski fc Kaluznvl 
ll99lT) so we applied straightforward aperture photometry 
using a large aperture of 3.5 arcsec which adds the flux from 
both stars and which was subsequently corrected (see next 
section). Several comparison stars within the field of view 
were checked for variability during each night and during 
the entire data set. 



3 THE R BAND LIGHT CURVES 
3.1 The ellipsoidal modulation 

The true V404 Cyg i?-band lightcurves were obtained after 
subtracting the contribution of the contaminating compan- 
ion star, assuming t hat it had i?=17.52±0.01 as given in 
ICasares' et alJ dl993l). We then folded our data using the 
ephemeris of ICasares fc Charles! lll994h . These lightcurves 
show the classical s econdary's ellipsoidal modulation (see 
Sha hbaz et al.lFl994h although with superposed variability. 

The ellipsoidal modulation is due to the differing as- 
pects that the tidally distorted secondary star presents to 
the observer through the orbital cycle. An ideal pure el- 
lipsoidal lightcurve is hence a double humped modulation 
with two equal maxima and two somewhat different minima 
(the minimum at phase 0.5 is deeper because gravity dark- 
ening is more pronoun ced in the direction of the L\ point 
jAvni fc Bahcalll975l) . This modulation was first simulated 
employing two sine waves at the orbital frequency and its 
first harmonic, where the fundamental allows the minima to 
vary. Relative phases were fixed to produce equal maxima. 
However, it has been largely reported that lightcurves of qui- 
escent tra nsients of t en ex h ibit unequal maxima as well (eg . 
A0620-00: iHaswelll Il996t iGelino. Harrison fc Oroszl l200lt- 
J1118 +480: IZurita et aTT 120021 - J2123-058: IShahbaz et alJ 
2003). The distortions from the theoretical case could be 
explained assuming a non-constant contribution of the ac- 
cretion disc, X-ray heating, spots on the surface of the sec- 
ondary star or light from the stream impact point. In any 
case, we also account for this possibility by allowing the sine 
phases to vary independently. We then fitted these two mod- 
els to the mean light curve and also to the lower envelope. 
The lower en velope was cons t ructed as has been previously 
explained in IPavlenko et alJ lll996l) and ZCS, that is, for 
each night we have selected all brightness estimates which 
differ from the minimum by no more than twice the accuracy 
of the observations. The selected data were then averaged 
in time and phase intervals. We found the best representa- 
tion (the minimum \ 2 ) fitting the lower envelope with the 
second model. These fits for the years with a most com- 
plete database (1992, 1998 and 2001) are shown in Figure 
Both minimum and maximum magnitudes varies. The max- 
imum at phase 0.25 is higher than at phase 0.75 and also 
the minimun at phase 0.5 is deeper than at phase 0. How- 
ever their relative strength changes from year to year. Apart 
from this, the peak to peak amplitude of the first harmonic 
remains roughly constant at 0.24±0.01 magnitudes, a key 
fact in the modeling of V404 Cyg and the determination of 
its compo nent masses. We note that the 1998 lightcurve ob- 
tained bv IPavlenko et al] i200 lT) in the same band exhibits 
equal maxima, however we cannot conclude whether this 
behaviour reflects a real change in the system, since both 
lightcurves are not simultaneous, or is due to an incomplete 
sampling. 



3.2 The superimposed fast variability 

We assumed that the fast variability superimposed on the 
ellipsoidal curve is not related to the ellipsoidal modulation 
so the former was isolated by subtracting the lower envelope 
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Table 1. Log of observations 



Night Band Resolution (s) Monitoring time (hr) 



1992 

lm JKT 



Jun 27 


R 


259 


1.45 


Jun 28 


R 


259 


5.44 


Jun 29 


R 


259 


2.51 


Jun 30 


R 


259 


7.11 


Jul 1 


R 


259 


7.07 


Jul 2 


R 


259 


5.33 


Jul 3 


R 


259 


6.18 


Jul 4 


R 


259 


5.15 


Jul 5 


R 


259 


4.11 


Jul 6 


R 


259 


6.16 


Jul 7 


R 


259 


6.37 


Jul 8 


R 


259 


6.10 


Jul 9 


R 


259 


7.24 


Jul 10 


R 


259 


5.67 


Jul 11 


R 


259 


6.68 


Jul 12 


R 


259 


6.29 


1998 








0.8 m IAC80/lm OGS 






17 Aug 


R 


94 


3.98 


18 Aug 


R 


88 


6.09 


19 Aug 


R 


79 


6.32 


21 Aug 


R/I 


129/99 


6.52/7.38 


22 Aug 


R/I 


146 / 99 


6.21/5.88 


23 Aug 


R/I 


136/129 


6.77/1.08 


24 Aug 


R/I 


136/129 


7.03/7.28 


25 Aug 


R/I 


138/129 


6.74/7.21 


26 Aug 


R 


136 


4.54 


27 Aug 


R 


136 


7.19 


1999 








lm JKT 








Jul 6 


R 


80 


5.89 


Jul 7 


R 


80 


6.05 


2001 








0.8m IAC80/lm TCS 






31 Jul 01 


R 


314 


3.80 


1 Aug 


R 


314 


4.12 


2 Aug 


R 


314 


5.22 


3 Aug 


R/J/K a 


314/1062/1100 


4.98/4.80/2.56 


4 Aug 


R 


314 


4.17 


6 Aug 


R/J/K a 


314/1054/1100 


3.88/1.80/3.72 


7 Aug 


R/J/K s 


314/1030/1040 


3.37/3.05/2.73 


8 Aug 


R/J/K s 


314 / 1030 / 1050 


4.04/2.53/3.83 


10 Augl 


R/J/K s 


314 / 1630 / 1050 


4.34/4.44/3.53 


11 Augl 


R 


197 


3.87 


28 Augl 


R/J/K a 


314/1100/1050 


5.19/5.32/5.05 


30 Augl 


R/J/K s 


314/1150/1050 


5.94/5.36/5.04 


31 Augl 


R 


314 


6.54 


6 Sep 


R 


314 


3.84 


8 Sep 


R 


314 


4.78 



variable-phase sinusoidal fit (as has been explained in the 
previous subsection) from the original data resulting in a 
light curve essentially free from the ellipsoidal modulation. 
In the case of the 1999 lightcurve, because only two nights 
of observations were available, we fixed all the parameters 
of the double sine wave obtained from the 1998 fit except 
the phase of the first sine. In 2003 only five nights of 
observations are available so we fixed the mean magnitude 
in the variable- phase fit at the value we obtained with 



the simplest constant relative phase fit. Furthermore, 
in order to compare the results from different years, we 
have resampled the lightcurves to have the same temporal 
resolution (320 s). The lightcurves and the corresponding 
lower envelope fits are shown in Figure |5] and the result of 
the subtraction is represented in Figure [3] 

Examining these lightcurves we find that the observed 
variations resemble the lightcurves of stellar flares, i.e. rapid 
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Table 1. Log of observations (continuation) 



Night Band 


Resolution (s) 


Monitoring time (hr) 


2002: 








lm OGS 








27 May 02 


R 


315 


Am 


zo iviay uz 




315 


3.90 


02 Jun 02 


R 


315 


4.27 


03 Jun 02 


R 


315 


4.42 


0.8 m IAC80 








01 Jun 02 


R 


315 


3.37 


04 jun 02 


R 


315 


2.05 


2003: 








lm OGS 








17 Jul 03 


R 


252 


1.20 


18 Jul 03 


R 


252 


8.20 


19 Jul 03 


R 


252 


4.33 


20 Jul 03 


R 


252 


8.26 


0.8 m IAC80 








28 Jul 03 


R 


352 


8.72 



and strong increases in brightness. We consider a flare to 
be any variation which differs from the mean level by more 
than 3cr. We found in many nights large single events of 
~2h or more duration appearing together with faster and 
small amplitude variations. Other nights, however, only the 
faster variations are present. The shape of the large flares 
seem to be symmetric with similar rise and decay times. In 
order to characterize the flares and quantify the flare activity 
we calculate the mean, z, standard deviation corrected for 
the instrumental error, a z (i.e. we subtract the instrumental 
error in quadrature ) and equivalent duration (fl z ) of the 
relative intensity lightcurves, i.e. relative to the steady non- 
variable level (see ZCS). 

Here z is the relative intensity, defined as: 

z = (I- / q )/J q (1) 

where I and 7 q are the observed intensities of the flaring and 
not flaring state respectively. 

The equivalent duration (II) is defined as: 




The total energy of the flare (E) is obtained by multi- 
plying the equivalent duration by the system's luminosity, 
i.e. E — 4 7r d 2 F,j II erg (where d is the distance to the 
system and F q is the dereddened observed quiescent flux). 
The time-averaged flare luminosity < Lr > is defined as 
the sum of the energies of the individual flares divided by 
the total monitoring time (Aft). The system properties and 
the resulting parameters of the variability are summarized 
in Tables |21 and respectively. 

We note that the flare mean level and the time-averaged 
flare luminosity are somewhat sensitive to the ellipsoidal 
model which is subtracted. It may be that the model is 
not an ideal representation of the non-variable level due to 
our incomplete sampling and effects other than the flares 
which distort the ellipsoidal modulation. In contrast, the 
standard deviation (a z ) is not significantly affected by the 
model assumed. We hence consider that a z is the best flare 



level activity indicator. We do not find any clear trend 
towards increasing or diminishing flare activity with time 
(see Table |H] and Fig. . Moreover, the long-term photo- 
metric study of V404 Cyg during quiescence performed by 
IPavlenko et alJ i200ll) also shows the fast variability per- 
sisting over the decade 1992-2000 with variable amplitude. 
In particular in 2000 the mean amplitude was 0.06. This is 
consistent with our data as is clear in Fig. 21 We estimate 
the flare activity to be roughly constant with z ~0.065, 
a z ~0.036 and < L R >~13x 10 32 ergs -1 . Assuming that 
the X-ray (0.3-7 kev) luminosit y is < L x >~5xl0 33 ergs" 1 
(from the Chandra spectrum; iKong et al.l 120021) we esti- 
mated the ratio between optical flare and X-rays luminosi- 
ties: < L pt/Lx >~0.3 



4 SEARCHING FOR PERIODICITIES 

As the short term variations are non-sinusoidal we have 
attempted to analyze and separate the periodicities present 
in the data by fitting a Fo urier series to the lightcu rves over 
a range of frequencies, v iMartinez fc Koenlll994l) . A plot 
of the reduced \ 2 °f each fit versus v should have a mini- 
mum at the fundamental frequency of oscillation in the data. 

We applied this method after resampling the data 
from different years to the same time resolution, and after 
detrending the long-term variations by subtracting the 
lower envelope as explained in the previous section. We 
rejected those nights having less than 3hrs coverage so 
as to avoid peculiar offsets due to poor sampling. The \ 2 
spectrum is computed in the range 3 to 10 cycles d _1 with 
a resolution of 5xl0 -2 cycles d _1 and with 3 terms in the 
Fourier series. They are shown in Fig. |5] where we have also 
marked the 68 percent confidence limit at Xmm + 1 ( e -g-> 
lLampton. Margon fc Bowvert 1976) 

With the 1992 data, we r eprod uce the results already 
published in IPavlenko et al.l il996J) . There exist a peak 
in the vicinity of 6hr (0 d .25) which is split into two (at 
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Table 2. V404 Cyg properties. 



Distance (kpc) 3.5±0.2 

A v 4.0 

Spectral Type KOiv 

inclination 56° ±2 

P orb (days) 6.64714 

Dereddened magnitudes R=13.5, 7=13.81, J=12.55, ii"=12.05 



Par ameters extracted frq m lWaemer et alJ lll992MShahbaz et alJ <1994MPavlenko et all il99rt : lMar"tinl il99fll) 
and ICasares et alJ il993l) . 



Table 3. V404 Cyg flares properties for the different years. 



Year 


< nt > 


Monitoring 


< L R > 


z f 


0> 






time (hr) 


xltfi 2 ergs -1 






1992 


0.055 


88.9 


10.7±0.4 


0.057 


0.040 


1998 


0.084 


61.4 


15.8±0.5 


0.084 


0.034 


1999 


0.066 


11.9 


12.5±0.4 


0.067 


0.025 


2001 


0.065 


68.1 


12.7±0.4 


0.061 


0.042 


2002 


0.071 


22.0 


15.1±0.5 


0.064 


0.038 


2003 


0.073 


30.7 


13.8±0.5 


0.067 


0.041 



(f) Time averaged equivalent duration of the individual flares divided by the total monitoring time. 



d .260 and d .246) with a central peak coinciding with 
the median value of these periods (at d .253). The lday 
aliases are also visible. A deta i led a nalysis of these data 
can be found in lPavlenko et ail <ll996T) . The 1998 spectrum 
shows two peaks at d .247 and d .233 exceeding the 68 
percent confidence interval and also their 1 day aliases. 
The power spectrum of the 2001 data is very noisy and we 
found two possible periods at d .241 and d .235 (also with 
their aliases) and two other peaks at d .261 and d .250 
all of them with the same significance. We note that these 
data were also analysed in ZCS. Here a Lomb-Scargle 
periodogram was performed yielding significant peaks only 
above 1 hour but with a complicated alias pattern which 
makes it difficult to distingui s h any periodicity. Finally, 
we also note that iHvnes et alJ <|2002T) analysis of the 1999 
lightcuves obtained at CrAO and SAI 1 together with 
our 1999 JKT data could not find any clear periodicity. 
Here the peak with highest significance has a period of 
d .254 although other periods are also possible. Clearly 
there is no coherent period in our database but variabil- 
ity on a characteristic timescale of about 6 hours is detected. 

5 THE COLOUR OF THE FLARES 

5.1 The 1998 simultaneous R and / photometry 

V404 Cyg was also observed simultaneously in the R and 
I bands during 1998. To extract the true I magnitudes of 
V404 Cyg, we previously subtracted the contribution of 
the contaminating star, I cs . We calibrated the unresolved 
pair of V404 Cyg plus contaminating star using the I 

1 with the 0.38-m telescope of the Crimean Astrophysical Obser- 
vatory (CrAO) and the 1.25- and 0.6-m telescopes of the Stern- 
berg Astronomical Institute (SAI) in 1999. 



magnitudes tabulated in lMartinl <ll999T) and then subtracted 
the real median magnitude of V404 to find I cs — 17.0±0.1. 
The resulting lightcurve is shown in figure |5] where fast 
variability is clearly superimposed on the ellipsoidal modu- 
lation. We detrended the long-term ellipsoidal modulation 
by subtracting the lower envelope as explained in section |3] 
We also binned the lightcurves to 300 s. 

Variability is clearly seen in the I band as well, 
following the same pattern as the variability in the R band 
(see Fig. We calculate the mean, standard deviation of 
the lightcurves and the equivalent duration of the flares 
relative to the steady non-variable level in both the R and 
I bands. We note that although the time resolution is the 
same, the parameters for the R band change from those 
reported in section HJ since now we are only including in our 
analysis the nights with silmultaneous observations. The 
results are given in table 2] 

The flare activity in the R band is larger than in the 
I band although we note that the relative contribution of 
the secondary star and the accretion disc to the total light 
is different at both wavelengths. In the case of V404 Cyg 
the distribution of the veiling rises toward the blue follow- 
ing a power law of index a — — 3.4 ± 0.9 which basically 
reflects the drop in the photosph eric flux of the last -type 
companion short ward of 5300 A iCasares et al.1 11993) . We 
also calculated the fluxes of the flares in both bands. To 
do it we first de-redden the observed magnitudes using for 
the reddening a value of Av=2.8 (see lShahbaz et al"l l2003b'l 
and then subtracted the ellipsoidal modulation constructed 
as explained in the previous sections from the de-reddened 
lightcurve yielding the monochromatic flux of the flares. We 
obtained Fr = 1.25 mjy and Fi = OAlmJy, so the mean 
flux density ratio of the flares is Fr/Fi ~ 3. Nevertheless 
there exists several sources of error in the calculation of the 
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Table 4. V404 Cyg flares properties for the R and / bands in 1998. 



Band 


< nt > 


Monitoring 


< L > 


Z S 








time (hr) 


x l(fi 2 ergs -1 




R 


0.097 


28.8 


18.2±0.6 


0.095 


0.030 


1 


0.038 


28.8 


3.8±0.1 


0.039 


0.019 



(f) Time averaged equivalent duration of the individual flares divided by the 
total monitoring time. 



flare fluxes. First, there is the photometric errors as well as 
the correction of the flux from the nearby star which contam- 
inates V404 Cyg (see section^. The photometric accuracy is 
~1 percent in both bands but the magnitude of the contam- 
inating star is known with a ~1% and ~10% accuracy in R 
and / respectively although, in fact, we must consider these 
values as lower limits since they were obtained from the lit- 
erature and there is difference between effective band passes 
on different telescopes. But the main source of error is the 
determination of the reddening. If we us e the extreme values 
for the reddening Ay =4.0 and Ay=2.2 llCasares et al Jl993t 
IShahbaz et alll 19941) we found an 80% error in the determi- 
nation of the ratio Fr/Fi. In addition we should note that 
the lower envelope could not reproduce accurately the sec- 
ondary's star lightcurve since there are several factors which 
distort the ellipsoidal modulation that we cannot quantify. 

To test if the flaring activity is simultaneous in the 
two bands, or if one lags the other, we calculated cross- 
correlation functions between them. This analysis identifies 
correlations and reveals the mean lag between the R and 
I varia bility. We used the Disc rete Correlation Function or 
DCF dEdelson fc Krolildll988^ . Our DCF, for the nights 
when simultaneous observations last more than 6 hours (21, 
24 and 25 Aug.), is shown in figure [7] Here positive lags 
correspond to R variations lagging behind those in /. The 
correlation is clearly significant at the 3a level although 
the peak appears asymmetric. This could correspond to 
R variations la s ting l onger than those in I. According to 
iTonrv fc Davis! (1979|), we estimate a delay of R with re- 
spect to / of (49±40s). The measured delay is clearly not 
significant at the 3 a level and is below the time resolution 
used, so is almost certainly not real. We hence conclude that 
both the R and / emission are simultaneous within the ob- 
servational errors. In fact this is the result we expect since 
both bands are close in wavelength and so probably originate 
in the same region. 

5.2 The 2001 simultaneous R and infrared 
photometry 

Because the fractional quiescent emission rises toward the 
infrared, it is of considerable interest to investigate the 
properties of the flare emission in this band. Even in the 
infrared the remaining contamination due to flares can 
have significant impact on the binary para meters derived 
from modeling the ellipsoid al variations (e.g. ISanwal et alJ 
119951: IShahbaz et all Il996ft . We also observed V404 Cyg 
simultaneously at optical and infrared wavelengths during 
2001. After performing aperture photometry around the 
target and its close co mpanion, we subtrac ted the latter 
by assuming J=14.53 llCasares et al] Il993f) and _K = 14.3 



iShahbaz et al] , 19961. The lightcurve is shown in figure |5J 
where fast variability is superimposed on the ellipsoidal 
modulation. 

We removed the long-term ellipsoidal modulation 
by again subtracting the lower envelope as in section |3] 
The non-ellipsoidal variability is shown in figure 1111 To 
obtain the same temporal resolution in all filters, we first 
binned the R data and then performed linear interpolation 
of the R and K s magnitudes to the times when the 
J magnitudes are given. We note variability above the 
noise level, also in the infrared bands. In some nights the 
infrared emission is clearly correlated with the optical (7, 
8, 10 and 28 Aug.), with amplitudes A J ~0. 16 mags and 
AK S ~0. 08 mags. On the night of 3 Aug. the K s flux 
shows a steep variation with 0.1 mags amplitude and ~1 
hour duration not correlated with the optical or J band flux. 

In this case we also obtained the flux densities as was 
made in the prev i ous se ction using the extinction law by 
iRieke fe Lebofskvl Q1985) . These fluxes in the regions with 
simultaneous photometry are fij=0.97mjy, Fj=0.68mJy 
and Fk =0.23mJy. Assuming Fr/Fi ~ 3 and a ~80 percent 
uncertainty in this ratio (see previous section) we estimated 
the flux density in the I band. Also we compare d our values 
with those obtained by IShahbaz et al. (2003b) to extract 
the spectrum of the flares (Fig 1101 . These last values are 
actually the fluxes at the peak of individual flares so we 
consider them as upper limits. Apart from this, we can only 
obtain a very rough estimation since we are comparing data 
from different epochs and with large uncertainties. 

In order to characterize the flare activity, we calculate 
the parameters of the variability in the three bands, the 
results being given in table We found lower activity 
at longer wavelengths although, as we will discuss in the 
next section, the relative contribution of the accretion disc 
depends on the waveband and hence a further correction 
must to be made. 



6 DISCUSSION 

Our R band lightcurves, obtained during the period 1992- 
2003, show considerable suborbital variability in addition to 
the ellipsoidal modulation. Assuming that these variations 
are not connected, we obtained a representation of the 
ellipsoidal modulation of the tidally distorted secondary 
star by fitting a double sine to the lower envelope of the 
lightcurve. The fits performed with the most complete 
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Table 5. Flares properties for optical and infrared in 2001. 



Band 


< nt > 


Monitoring 


< L > 


*/ 








time (hr) 


xltP 2 ergs -1 






R 


0.071 


21.1 


13.5±0.4 


0.069 


0.036 


J 


0.057 


21.1 


7.7±0.3 


0.058 


0.029 


K 


0.035 


21.1 


1.8±0.1 


0.029 


0.012 



(f) Time averaged equivalent duration of the individual flares divided by the total monitoring time. 



databases (years 1992, 1998 and 2001) show unequal 
maxima and minima with relative strength varying from 
year to year. This behavior could be explained assuming 
starspots on the secondary star. This is a very widespread 
phenomenon among cool stars. Axial rotation of a spotted 
star and slow variations of the starspot geometry cause 
photometric variability which is expressed in brightness 
rotational modulation and slow variations of mean stellar 
lightcurves. Other possibility is they are produced by 
changes in the accretion disc geometry and brightness due a 
superhump, i.e. optical modulation as a consequence of the 
precession of an eccentric accretion disc by perturbation 
from the secondary. This is a very promising suggestion 
since superhumps in a quiescent SXT have already been 
detected fJ1118+480. IZurita et"ai]|2002ft . Superhumps have 
also been largely invoked to explain the changing quiescent 
lightcu rve of the prot otypical black hole candidate: A0620- 
00 fee. lHa swcll 199^). After subtracting the corresponding 
lower evelope fits from the original data the lightcurves are 
dominated by large flares (~ m .2) lasting a few hours. In 
addition, shorter timescale variability (~ m .05) is present 
with lower amplitude. Our photometry indicates that the 
flaring activity did not change significantly over this 10-year 
interval which indicates that this variability is probably not 
connected to the 1989 outburst. 

We have also found significant variability in the I and 
near-IR band (J and K s ) and these are correlated. In 1990, 
just one year after the outbu rst, a variation of ~ m .2 was 
seen in several visible bands (iCasares et alJll99St) . In 1991, 
the variability originally interpreted as a six hour modula- 
tion, decreased to ~ m .l, but no significant cha nges were 
reported in the interval 1992-1993 <MartirJll99^ . Infrared 
photometry in 1990, was dominat ed by a ~ m .2 amplitude 
at a periodicity of ~6 hours (ICasares et aljll993l) . Also 1993 
H band data presents signific ant variability contam inating 
the ellipsoidal modulation dSanwal et alJ Il995l) . They 
reported an rms scatter of ~ m .04 for V404 Cyg and 
~ m .02 for the comparison star. This yields a flare activity 
level of <J Z -~0.032, which is consistent wi th our values. 
In 199 3, K band photometry obtained by IShahbaz et alJ 
il994f) does not show short term variability, with an upper 
limit of ~ m .03. However in this case, only one night of 
more than 6 hours of continuous coverage was obtained. 

We have also found that lower activity is found at longer 
wavelengths. However, as we have pointed out in ZCS the 
flare activity level calculated as <r z represents the flare activ- 
ity relative to the underlying smooth modulation (F q ). F q 
contains two distinct components: the flux from the com- 
panion star, F c and the flux from the non-variable, steady 



accretion disc, F<j (see figure 9 in ZCS for clarity). The rel- 
ative contribution of the accretion disc depends on the wave 
band and hence, if we ass ume that the flare s are produced 
in the accretion disc fsee lHvnes et alJl2002l . ZCS), a more 
useful parameter to compute is the flare activity relative to 
the flux from the steady disc, rather than relative to the 
total non-variable flux, which is given by 



ffz = o"z/r?d 
where 

Va = F d /F q = (l + F c /F d y 1 



(3) 



(4) 



We can determine rj d using the spectroscopic value for 
the veiling ?7obs, which is defined as 



Vobs 



F f + F d 



(•») 



Ft + F q 

where Ff is the mean flux of the flares, and find 

rjd = (zf + 1) ?7obs - Zf (6) 

ICasares et alJ lll99St) reported that the accretion disc 
contribute s 11-16 percent of the observed flux in R (see 
also ICasares fc Charles! <ll994h ). Therefore we assume an 
observational veiling of r^obs = 13% in this band. We 
also estimate an observational veiling of ~ 12% in the I 
band by extrapolati ng from the veiling in B, V and R 
( Oasares et al.lll995t . On the other hand, IShahbaz et alJ 
determined the fraction of the light arising from the 
accretion disc in the K band varies from to an upper limit 
of 11 percent. We hence consider 77obs=5 percent in this 
case. The corrected values using these factors, are tabulated 
in table HJ The uncertainties on the veiling and the limited 
database only allows us to say that the activity diminishes 
slowly or remains approximately constant through the IR. 
In the 1992 series of B, V, R and / band observations (16 
measurements spread over 16 nights; see lMartirjH99"gh an 
estimate of the variability could be obtained by comparing 
the variance with the photometric noise and the optical 
variations. This variance increases toward the blue suggest- 
ing that the source of fa st variability in V404 is indeed blue. 
We note however that IShahbaz et alJ (I2003IJ) determined 
the colour of the large flares and found the i' (A e //7650A) 
band flux to be larger than that in the g' (A e //4750A) 
band. Comparing with these data we obtained an estimated 
spectrum with a peak at A e //7000A(see Fig. 1101 . 

In an attempt to compare the infrared variability 
we also observed the SXT with highest flares: J0422+32, 
where events as high as 0.6 mag where observed (ZCS). 
Observations of this system in the H band covering half of 
its orbital period do not show variability above ~0. 06 mags, 
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Table 6. Flare activity and veiling factors for R, I and K band. 



Band 


Tfoba (percent) 


rjd (percent) 


*/ 








R 


13 


14.1 


0.095 


0.030 


2.1x10" 


-3 


I 


12 


12.4 


0.039 


0.019 


1.5x10" 


-3 


R 


13 


13.8 


0.069 


0.036 


2.6x10" 


-3 


K 


5 


5.1 


0.029 


0.012 


2.3x10" 


-3 



the estimated error of our photometry (Zurita et al. in 
preparation). Comparing with the variability we observed in 
the R' band (ZCS) and assuming that the disc contributes 
~30% in H (Beekman et al. 1997) we estimate that the 
activity in the infrared is, at least, 33% lower than in the 
optical. Although only a very crude estimation can be 
made, among other things because the observations are not 
simultaneous and because the value of the veiling is not 
well determined, it seems that the spectrum of J0422+32 
flares also drop at longer wavelengths. On the other hand, 
no significant contamination of the infrared lightcurv es of 
J0422+32 has been found bv lGelino fc Harrison! (12003ft . 

The quasi pe riodic oscillations not iced in the earliest 
data (1992 - see IPavlenko et al.l Il99fit) are less clear in 
our more recent data (1998 and 2001). However, the 
periodograms show peaks close to 6 hr at the 68 percent sig- 
nificance level which could be interpreted as a consequence 
of the appearance of a flare event almost every night. 
Shah baz et alJ l)2003Fl found a QPO feature at 21.5 minutes 
attributed to variability at the ADAF transition radius. 
Our longer sampling does not allow us to test this hypoth- 
esis. It is possible in any case that variability with different 
timescales have different origin. 
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Figure 1. Ellipsoidal model fits to the lower envelopes of the 1992 
(continuum line), 1998 (dashed line) and 2001 (pointed line) R- 
band lightcurves. 
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Figure 2. The R band ellipsoidal lightcurves of V404 Cyg for 
1992, 1998, 1999, 2001, 2002 and 2003 and the lower envelopes 
constructed as explained in section|!J] Different magnitude offsets 
were aplied for clarity. 
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Figure 3. R band lightcurves of V404 Cyg after subtracting the 
ellipsoidal modulation - 1992 lightcurves. 
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Figure 3. i? band lightcurves of V404 Cyg after subtracting the 
ellipsoidal modulation - 2001 lightcurves. 
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Figure 3. i? band lightcurves of V404 Cyg after subtracting the 
ellipsoidal modulation - 1998 lightcuurvcs. 
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Figure 3. R band lightcurves of V404 Cyg after subtracting the 
ellipsoidal modulation - 1999 lightcurves. 



CD 



0.3 

0.2 

0.1 


0.3 

0.2 

0.1 





I I I I I I I I I I I 

27/5 J 

->""• •' 

.* 

\ 3 cjj _ 


1 1 1 1 1 1 1 1 1 1 1 

28/5 - 

. ' i . 




I | I | I | I | I | I | I 
1/6 - 

vVT; 3 

i " 


i i i i i 

2/6 - 

X • * 3 up 




I i I i I i I i I i I 
3/6 

« • 




I i I i I i I i I i I 

4/6 - 

3,^ 





012345601234560123456012345 

Time (hours) 

Figure 3. Continued - 2002 lightcurves. 
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Figure 3. Continued - 2003 lightcurves. 
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Figure 6. The R and I 1998 lightcurves and the lower envelopes 
constructed as explained in section 131 
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Figure 9. The J and K s ellipsoidal 2001 lightcurves and the 
lower envelopes constructed as explained in section 131 



2 pr 

x 0.8 — 

£ 0.4 - 

L 



g Ri I J 



k- 



0.5 1 1.5 

X(p, m) 



2 2.5 
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Figure 8. The 1998 R, I and color R-I lightcurves of V404 Cyg 
after subtracting the ellipsoidal modulation. 
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Figure 11. From top to bottom: The R, J and K s lightcurves of V404 Cyg (after subtracting the ellipsoidal modulation) and the color 
indices R — K s , J — K s and R — J. Crosses mark the original data and circles the new resamplcd data. 
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